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In this paper we examine thoroughly Higgs to /xr+r/x decay via processes involving 
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R-parity violating parameters, the branching ratio could still achieve notable value 
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interest to future experiments. We present here key results of our analysis. Based 
on the analysis, we give some comments on h° — > ep, + /ie and h — > ef + re also. 
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I. INTRODUCTION 



As we know, in the Standard Model (SM) lepton number of each flavor is separately 
conserved. Thus lepton flavor violating (LFV) decays such as Higgs to /if + t/J is forbidden. 
However, neutrino oscillation experiments provide strong evidence that the lepton flavor 
conservation should be violated |lH4|. If lepton flavor violation can be observed in processes 
involving only SM particles, this would contribute an important probe to physics beyond 
the SM. Such processes, in particular the Higgs to /if + rp, decay, deserve attention. 

Looking into the literature, various sources or scenarios to accommodate LFV interactions 
have been introduced and analyzed. For examples, adding heavy right-handed neutrinos 
can give neutrinos mixings and hence lepton flavor violation j^|, or a general two Higgs 
doublet model has LFV interactions due to Yukawa coupling matrices which could not be 
diagonalized simultaneously jfilQ]. Under the framework of supersymmetry (SUSY), it is well 
known that nonzero off-diagonal elements of soft SUSY breaking terms in the leptonic sector 
(A , fh 2 L and fh 2 E - to be precisely defined below) generate LFV couplings. Moreover, the 
(total) lepton number itself may not conserve. For the SUSY case, such R-parity violating 
(RPV) couplings also give interesting contributions to processes like Higgs to /if + r/2. 

While SUSY is undoubtedly a popular candidate theory for new physics, its existence is so 
far in lack of experimental evidence 8] . Thus, some simple versions of supersymmetric model 
such as constrained minimal supersymmetric standard model, have faced stringent challenges 
{9]. However, it is pointed out that there are still rooms for (minimal) supersymmetric 
standard model to accommodate existing experimental constraints 10Ml3|. For instance, 
large mass spectrum for the majority of supersymmetric particles around or beyond ITeV 



is still to be probed 
boson mass = 125GeV to 126GeV 



13j. The heavy spectrum is in accordance with the newly discovered 



14 



16] . A large portion of the parameter space remains 



uncovered in versions of the minimal supersymmetric standard model (MSSM) with more 
free parameters Non- universality of soft SUSY breaking masses is also a possible 

explanation to the non-observation of supersymmetric signals 12]. 

Under the scheme of the MSSM, various LFV decays such as r — > /I'j, r — > fiX, r — > fir], 
t — > fififi, and so on [l^], as well as Higgs to fif + rp, decay [3, Q] which we put our focus on 
in this paper, have been discussed. However, in many studies of the MSSM, R parity is often 
imposed by hand to prevent proton decay and make the lightest supersymmetric particle a 
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possible dark matter candidate. From the theoretical point of view, R parity is ad hoc and not 
well motivated so long as the phenomenological (minimal) supersymmetric standard model 
is concerned 19]. A generic supersymmetric standard model (without R parity imposed), 
on the contrary, not only provides a convenient way to lepton flavor violation, but also has 
the advantage of a richer phenomenology including neutrino masses and mixings without 



introducing any extra superfie 
there have been some studies 



d. Under the framework of SUSY with R parity violation, 



20 



21] on the issue of lepton flavor violation. Nevertheless, 



such studies were either limited to particular types of R-parity violation or did not take 
h° — > fir + T/i into consideration. While recently both ATLAS and CMS [jj], [l5] reported 
discovery of a boson state which is essentially compatible with a SM-like Higgs, more data 
is needed to pin down its nature and the flavor violating Higgs decay such as h° — > \if + rft 
is especially interesting at this moment. In this paper, we will investigate thoroughly LFV 
Higgs to jif + tJi decay from SUSY without R parity via full diagrammatic calculations 
up to one-loop level. Best branching ratios for the decay from various RPV parameter 
combinations, under a reasonable choice of an experimentally viable region of the background 
parameter space, will be reported. Note that a part of key results has been reported, with 



22]. 



limited presentation of analytical expressions and discussions, in a short letter 

In following section, we summarize our basic formulation and parameterization of the 
generic supersymmetric standard model (without R parity). Particularly, the neutral and 
charged Higgs mass terms, including loop corrections (up to two-loop for neutral Higgs case), 
would be discussed. Then we give a sketch of our calculations and show numerical results 
from all possible RPV parameter combinations in section 3. Note that during our analysis 
we made no assumptions on the RPV parameters. The mass spectrum of all SUSY particles 
as well as the Higgs boson are kept within experimental constraints. Finally, we conclude 
this paper with some remarks in section 4. Lists of all one-loop diagrams and useful effective 
couplings would be given in appendices in the last. We may be including more details than 
necessary, particularly in the sense of showing some experimentally uninteresting results. We 
include those to give a full picture about the physics involved, so that readers can appreciate 
the key features leading to the interesting or uninteresting results. Some of the lessons one 
can learn from the analysis would be useful for future studies of other related aspects of the 
model. Under the same consideration, we give detailed expressions of the couplings involved 
and the Feynman diagrams in the appendices. 
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II. ON THE SUPERSYMMETRIC STANDARD MODEL WITHOUT R PARITY 

AND THE SCALAR MASS MATRICES 

A. Formulation and parameterization 

With content of minimal superfields spectrum, the most general renormalizable superpo- 
tential without R parity can be written as 



W = e ab 



where (a, b) are SU(2) indices with 612 = — €21 = 1, (i,j,k) are the usual family (flavor) 
indices, and (a, (3) are extended flavor indices going from to 3. Note that A is antisymmetric 
in the first two indices as required by SU(2) product rules while A" is antisymmetric in the 
last two indices by SU(3) C . The soft SUSY breaking terms can be written as follows: 



+ ^: jk U?D?D% (1) 



2 



+ h.c. 



Koft = e ab B a HILl + e ab A^H^ + A^ d Q\D) + A^H a d L\E 

1 

2 

+ Q ] rh 2 Q Q + U^rhfj U + D^m 2 D D + Dm 2 L L + &m\ E + m 2 Hu \H U 



+ Cafe 



+ l^OPDfb? + h.c. 



+ ^BB + fwW + fsS + ^, (2) 

where DmyL is given by a 4 x 4 matrix. m| oo corresponds to fn 2 Hd in MSSM, while m| ofc 's 
give new mass mixings. 

The above, together with the standard (gauged) kinetic terms, describe the full La- 
grangian of the model. We have four L superfields, which contain the components of fermion 
doublet as 1° and l~, while their scalar partners as 1° and /~. In principle, the neutral scalar 
part l° a of all four L superfields can bear vacuum expectation values (VEV). To make the 
analysis simple and physics more transparent, we use a parameterization which picks a basis 
such that the direction of VEV is singled out, i.e. only L bears a nonzero VEV among 
four L's. This procedure guarantee Lq can be always identified as Hd in MSSM. The two 
superfields have the same quantum number as the symmetry of lepton number which makes 
the distinction between L and Hd as by definition not part of the model. However, one 
should keep in mind that Hd may contain partly the charged lepton states. It is also worth 
mentioning here that the down quark and charged lepton Yukawa coupling matrix are both 
diagonal under our parameterization while up quark Yukawa coupling is the product of CKM 
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factors and diagonal quark masses. The parameterization has the advantage that tree level 
RPV contributions to fermion mass matrices are described completely by the /ij, and 



fh 2 Loi parameters, which are well constrained to be small even with just very conservative 
neutrino mass bounds imposed [23 



|24fl. 

Now we turn to the issue about mass matrices of matter fields. In our framwork, the three 
known charged leptons, together with two charginos, correspond to the mass eigenstates of 
a 5 x 5 charged fermion matrix Aic-, which can be diagonalized by two unitary matrices as 
V*M C U = diag {M x - } = diag {M cl , M c2 , m e ,m At ,m T }. For neutral fermions, we take four 
heavy neutralinos and three very light neutrinos as mass eigenstates under the scheme of a 
7x7 neutral fermion mass matrix Aij^. By using a unitary matrix X, the diagonalization 
can be done as X T M.j K fX = diag |M x o j = diag {M„ fel4 , m Ul , m V2 , m V3 \. On considering 
the squark sectors, the up-squark mass-squared matrix looks exactly the same as the one in 
MSSM, while down-squark one contains new contribution from RPV terms. They can be 
diagonalized separately as V U ^M 2 V V U = diag{A4|} and V d] M 2 D V d = diag{Ml}. A11 the 
mass matrices mentioned above can be found in 



23] 



B. Scalar mass matrices and loop corrections 

For neutral scalar mass matrix, we have now five neutral complex scalar fields from H u 
and four L a 's. Explicitly, we write the (1 + 4) complex fields in terms of their scalar and 
pseudoscalar parts, in the order {h$ , 1$, Z°, I®, l®} to form a full 10 x 10 (real and symmetric) 
mass-squared matrix, which (in tree-level) can be written as 

„ / Ml s Ml P \ 
M 2 S = ss SP , 3 

\{Ml P ) T M' PP ) 
where the scalar, pseudoscalar, and mixing parts are 

M 2 SS = Re(Ml) + 2M 2 ^, 
M 2 PP = Re(M 2 ) , 

Ml P = -Im(Ml) , (4) 
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respectively, with 



Mt 



-Ml 
2 z 



sin 2 /3 


— cos/3 sin/? 


1X3 


cos/3 sin/? 


cos 2 /3 


1X 3 


O3X1 


O3X1 


O3X3 



(5) 



Mt 



™?h u + »>« + M 2 Z cos2/3 [-1] 
-(^) 



(6) 



and 

-(£«) 

™l + to/0 + ^1 cos2/3 [i] /, 

As for charged (colorless) scalars, we should treat charged Higgs and sleptons on an equal 
footing. The basis {hp, Iq , , 1%, lp , lp JP} as 1 + 4 + 3 form from H u , four L Q 's and 
three Ef 's is used to write the 8x8 charged scalar mass-squared matrix, which can be 
written as 

/ A^2 A~7 2 t A/, 

l RH 



Ml 



1 ^ ^ H ^ ^ 
M\ H Ml L Mr L 

M\ H M\ L M\ R j 



(7) 



where 



M 2 H 



M 2 



LL 



M 2 RR 



™>h u + f£l*a + Ml cos2/3 



- sin 2 fl 



w 



+ Mf sin 2 /3 [1 - sin 2 # 



w\ 



+ 



I + + (vM + M l cos2/3 

Mf cos 2 /3 [1 - sin 2 ^] 1X3 



1 



+ sin 2 ^ 



3xl 



O3 



= m| + m E m} E + Mf cos2/3 [- sin 2 ^] J 3x3 , 



(8) 



and 



M\ u 



M\ H 



(M 2 RL ) T 



\ M 2 Z sin2/3 [1 - sin 2 ^] 

O3X1 



-(K^ofc) ^= = 







* \ \ u 



(no sum over k) , 



(9) 
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with m L = diag{0,mE} = diag{0, m 1; m 2 , m 3 }. m,;'s (~ m e . under small-/ij scenario) are 



mass parameters in charged fermion mass matrix [23J. Furthermore, the two scalar mass- 
squared matrices can be diagonalized as V sT M 2 s V s = diag{Mj m=1 10 } and V lf M 2 E V l = 
diaejM? |, which will become useful later. 

Different from MSSM, the physical scalar states are now mixture of Higgses and sleptons. 
The RPV terms provide new contribution to scalar mass matrices and hence Higgs mass. In 
addition, radiative corrections, especially those from third generation quarks and squarks, 
could play important roles in Higgs mass. Accordingly, we implement complete one-loop 
corrections [25] to matrix elements directly relating to Higgs bosons (CP-even, CP-odd 
and charged ones as well) during our computation. Moreover, light Higgs mass should 
be treated delicately because of the newly discovered boson mass = 125GeV to 126GeV by 



Large Hadron Collider (LHC) [14l.ll5|. Therefore we include further an estimation 26| of key 
two-loop corrections in light Higgs related elements 1 . Note that radiative RPV corrections 
are typically too small to be taken into account, thus we study tree level RPV effects only. 

Under the scheme of MSSM without R parity, the one-loop effective Higgs potential is 
(recall that H d = Lq after our parameterization being chosen): 

V eS = (m 2 Hu + \fi a \ 2 ) \H U \ 2 + (m 2 Lm + \no\ 2 ) \H d \ 2 + (e ab B H«H b d + h.c.) 

+ I (9l + 9' 2 ) \H U \* + I (gl + 9 ' 2 ) \H d f + \ (g 2 - g' 2 ) \H U \ 2 \H d \ 2 - \g\ \e ah Hluf 



3 >. ^ 

+ 32vr 2 ^ 

q=t,b 



<H I Q2 

i=l,2 V ^ 



•■' 2 3\ /, ml 3 



. - 2mt( In ^| - - 
2 J q \Q 2 2 



(10) 



where Q is the renormalization scale which should be around weak scale (10 2 to 10 3 GeV). 
m 2 and m 2 denote (Higgs background fields-dependent) eigenvalues of the squark and quark 
mass matrices respectively. 

By using following linear expansion of Higgs bosons (with a relative complex phase for 
generality) , 

H,= \ K ), H d = e* ( 7! ^ + h ' d + ih " d} I , 

J, (,,„ + K - i K ) j y h- 



1 Though the Higgs bosons mix with the sleptons via RPV terms, we can still identify the Higgs bosons 
out of other sleptons due to the foreseeable smallness of RPV parameters. 
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tadpole equations can be written as 



v d Re(B e 



3 



+ 



16vr 2 



E E m ! 

.g=t,6 i=l,2 



- 2m; 



dm 2 
~dht 



; u Re(£ Q e^) = (m£ 00 + |/i | 2 ) v d + - (<? 2 2 + </ 2 ) u d (t; 2 - t; 2 ) 
3 



16vr 2 



E E ™l 

,q=t,bi=l,2 



dm 2 n . \ l m"-, 



1 - 2mf 



9m 2 
0# 



In 



Q 2 



- 1 



i0>. 



+ H^EE'«l /a "' 



q=t,bi=l,2 



(12) 



where m~ = (m 2 



) is the squark mass-squared, while expressions for the derivatives with 



respect to the scalar fields in the bracket (inc 
plicated so we do not list here. One can see, 



uding second derivatives used later) are com- 



25] for example, for details 2 . 



Tadpole equations along the direction of other scalars/sleptons can be obtained easily 
from scalar potential terms which are related to neutral sleptons 3 : 

y=J2[ + tf^) + (™L> + *>ti) h °JT + K w + fit! 



while vanishing derivatives of V give 



\h°\ 2 + 2 



0|2 



Bihffi + h.c. 



(13) 



Bi tan /3 = m 2 w + //JJ/a 



(14) 



The exact form of tree level elements of scalar matrices are as mentioned above, while 
the one-loop corrections from third generation quarks and squarks are: 

3 



M 



Loop 



3k 167r 2 



EE 



q=t,b i=l,2 



dm 2 \ /dm 2 , \ m~. „ / d 2 m 2 
%l \ I 'h. \ In 'h. -urn 2 - ' 

0fc \d<j>J Q 2 + m * 



- 2 



dm 2 \ /dm 2 \ m 2 
dfi / \d<t> k / Q 2 "Xd^dfa 



d 2 m 2 



1 



2 



(15) 



2 There may be sign difference between the expression for derivatives in the reference and ours due to 
definition of linear expansion of scalars. 

3 Conceptually, If is not the usual Vi since If deviates from Vi slightly, with parameter fii characterizing the 
deviation between them. See [23| for details. 
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In the case of neutral scalars, j and k can be any number among 1,2,6,7 which is corre- 
sponding to h s u) h s d , h a a and h a d respectively. As to charged scalar case, j and k can only take 
the value of 1 or 2, with 4>j = h d } and 4>k = {h~, h d }. With the one-loop corrections 
and the estimation of two-loop corrections 26[ added to the mass terms, we can make sure 



the Higgs mass generated by our program being with enough accuracy. 

III. CALCULATIONS AND NUMERICAL RESULTS 

At tree level, neutral Higgs can decay into /if or rp directly via RPV neutral scalar- 
charged lepton-charged lepton coupling which is absent in MSSM. Otherwise, it can decay 
through one-loop diagrams (or higher loop diagrams). We list in appendix A all possible 
one-loop diagrams containing RPV couplings for a neutral scalar decaying to [if while RPV 
effective couplings we used among all relevant mass eigenstates are listed in appendix B. In 
our numerical computation, we put all the tree level mass matrices into program and deal 
directly with mass eigenstates which is obtained by diagonalizing corresponding mass matrix 
numerically. The one-loop and two-loop corrections as mentioned in section 2 to matrix 
elements which is most relevant to Higgs mass are also implemented. To make analysis 
representative enough, we have fully calculated all the (tree and one-loop) diagrams that 
may contribute. By encoding the analytical formulas of decay amplitude into program, and 



using LoopTools [27| program for the evaluation of loop functions, the numerical value of 
total amplitude and hence decay rate can be obtained. With total width being the RPV 
decay rate of h° — > pf + rp we got, plus MSSM one which is obtained by including all the 
widths of significant decay channel such as bb, t~t + , WW*, ZZ*, 77, and gg, we have the 
branching ratio of h° — >■ pf + rp. 

Our aim is to use a concrete setting that is compatible with known constraints but 
not otherwise too restrictive, to illustrate what we expect to be more generic features of 
the RPV signature. During numerical analysis, following restrictions and assumptions are 
used: |/i |, \A U \, \A d \ and |A A | < 2500GeV; A e is set to be zero since its influence is 
quite small and negligible; m 2 E = m 2 L (without zeroth component) < (2500GeV) 2 with 
off-diagonal elements zero; tan/3 ranges from 3 to 60; heavy Higgs/sneutrino and charged 
Higgs/slepton masses are demanded to be between 200GeV to 3TeV. Furthermore, after 
considering the uncertainties in experimental Higgs mass and loop corrections to Higgs 
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mass terms, we kept numerical light Higgs mass to be in the range of 123 to 127GeV. We 
adopt the relation M 2 = = 1M\ between three gaugino masses and the condition 

that squarks of the first two families cannot be lighter than about O.8M3 . Therefore we 
take soft SUSY breaking scalar masses rhq = m\j = fh 2 D = (O.8M3 x identity matrix) 2 for 
simplicity and M 2 < 2500GeV in our analysis. The parameter setting is in accordance with 
the gravity-mediated SUSY breaking picture 



28 



for instance. 



There are several different sources which can give constraints on our RPV parameter 
setting. Among which, the one from indirect evidence of neutrino mass, i.e. Y2i m vi ~ 
29I ] is quite crucial. However, since it has not been completely ruled out for neutrinos having 
mass larger than leV, we give some comments on branching ratios as references under the 
condition that neutrino mass constrained only by the solid bounds, i.e. m Ue < 3eV, m Ufi < 
190KeV and m UT < 18.2MeV [30] as well. Besides, to prevent Bi from being unreasonably 
large, we set by hand 1% of B as upper bound in the circumstance of extraordinary loose or 
no available experimental bounds. Note that all LFV couplings/mass mixings that conserve 
R parity have been turned off during our analysis. That is to single out the effects of the 
RPV ones. The reported numerical branching ratios are the best numbers we found under 
the framework. 



A. Contribution from Bi Bj combinations 

The constraints on this type of combinations are mainly from neutrino mass experiments. 



The RPV parameter Bi can give contributions to neutrino masses via one-loop diagrams |31 |. 
Generally speaking, larger sneutrino and neutralino masses will raise the upper bound of Bi. 

Except for the combinations B 2 B% and B1B1, BiBj combinations can only give contri- 
butions to the decay from the Type2No.4 diagram (in Appendix A). Since heavy charged 
scalar masses will severely suppress this diagram, we can have relatively larger amplitudes 
only in the existence of light charged scalar (s). For the B 2 B 3 combination many diagrams 
contribute, hence its behavior is quite complicated. Basically, the decay amplitude from 
B 2 B 3 tends to increase when soft SUSY breaking scalar masses and gaugino masses get 
heavier due to the rise of the upper bound on Bi from neutrino masses as mentioned above. 
Note that B 2 or B3 alone can give contributions to the h° — > \if + rfx decay as well. Such 
contributions are, unavoidably, included in all combinations containing B 2 or B%. Combina- 
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tions that are not listed give zero contribution at one-loop level. Similar for the other kinds 
of RPV parameter combinations given below. Our results are shown in Table 1. 



Table 1. Bi Bj contributions to branching 


ratio of h° — > [it + t/j, 


RPV Parameter 
Combinations 


With Neutrino Mass 
<leV Constraint 


With Relaxed 
Neutrino Mass Constraint 


B\ B 2 


4 x KT 22 


3 x KT 21 


B\ B 3 


3 x KT 22 


3 x KT 21 


B2 B2 


9 x KT 23 


4 x KT 15 


B2 B 3 


2 x KT 11 


3 x KT 10 


B3 B 3 


8 x KT 23 


4 x KT 16 



B. Contribution from Bi fij combinations 

The Bi fij type of combinations gets constrained from several sources. The values of Bi 
and Bijij are highly constrained separately by their loop contribution to neutrino masses 



3jJ . On the other hand, a non-zero fij will induce tree level neutrino mass, hence it is 
also constrained. The non-observation of leptonic radiative decays like fi — > e^y, etc. also 
give upper bounds on Bfjju, say, \Bl{i 3 \, | -S2 A*3 1 5 1-^3^*1 an d I-B3/X2 1 ~ 10~ 4 |/io| 3 ; \B*^2\ and 



1^*1 < 7x 10- 7 |^o| 3 [2l|. 

All Bifij combinations except -£>2/-t3 and -83^2 can give contribution to the Higgs decay only 
from the Type2No.4 diagram. Again light charged scalars preferred for the case. However, 
these contributions (from the Type2No.4 diagram) can not provide a significant branching 
ratio. Hence we have the uninterestingly tiny numbers as shown in Table 2. 

As for -82/^3 and -83/12, both giv e contributions to the Higgs decay from several diagrams. 
Among them, the tree diagram [22| is the most important one over a wide range of parameter 
space. Especially, for the -83/12 combination, a key contribution to the decay amplitude is 
enhanced by the tau Yukawa coupling y e3 via a term proportional to 

y e3 M*B 3f i* 2 / [(/i M 2 - M^sin2/3)M S 2 ] 

(M 2 denotes a generic real scalar mass eigenvalue). The later makes the branching ratio 
from -B 3 /i 2 the largest one among all -Bj/i^'s. There is a similar feature for the contributions 
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from the B 2 fi 3 combination, but with a muon Yukawa y e2 instead. These two combinations 
get their best values under small fiQ and M 2 as can be seen from the expression above. 

Table 2. Bi fij contributions to branching ratio of h° — > fif + Tp, 



RPV Parameter With Neutrino Mass With Relaxed 

Combinations ^$l e V Constraint Neutrino Mass Constraint 



Bi fi 2 


1 x KT 24 


6 x KT 22 


B\ fi 3 


1 x KT 24 


6 x KT 22 




9 x KT 23 


8 x KT 22 


B 2 iii 


4 x KT 26 


6 x KT 15 


B 2 il3 


1 x KT 15 


9 x KT 6 


B 3 fii 


8 x KT 23 


7 x KT 22 


B 3 iii 


1 x KT 13 


7 x 10- 4 


B 3 il3 


4 x KT 26 


4 x KT 16 



If we allow neutrino masses being larger than leV while making jj large enough to 
unfasten restrictions by leptonic radiative decays, the main constraint will be from B^fij loop 
contribution to neutrino masses. Specifically, by demanding B2H3 and -83^2 contribution 
to neutrino masses to be just smaller than 190keV, their corresponding branching ratios 
increase significantly to the order of 10~ 6 and 10 -4 . Note that as in BiBj case, contributions 
from B 2 and B% alone are included. 

In fact, analyses similar to the above can be applied to h° — > ejl + [ie and h° — > ef + re as 
well. The Bifij contributions to h — > eft + fie is expected to be tiny due to the smallness of 
the corresponding Yukawa couplings y ei and y e2 . On the other hand, while the contributions 
from Bifi 3 is also suppressed by a relative factor of y ei /y e2 , and the contributions from B 3 fii 
to h° — > ef + re could be roughly the same order as that of h° — > fif + rji. Hence, the 
h° — > ef + re decay may also be of interest. 



C. Contribution from B; A combinations 



Apart from the constraint on the Bi parameters, the A type parameters are bounded by 



charged current experiments 32J]. Generally speaking, increasing soft SUSY breaking slepton 
masses and gaugino masses leads to heavier chraged slepton, sneutrino and neutralino masses 
hence raises the upper bounds for Bi and A. 
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Table 3. Bi A contributions to branching ratio of h° — > [it + r/i 



RPV Parameter With Neutrino Mass With Relaxed 

Combinations ^$l e V Constraint Neutrino Mass Constraint 



B\ A123 


1 X 


IO -5 


4 X 


IO" 5 


B\ A132 


3 x 


IO -5 


7 x 


IO -5 


B\ A232 


4 x 


1Q -22 


1 X 


1(T 21 


B\ A233 


5 x 


1Q -25 


2 x 


10- 24 


-S2 A 123 


7 x 


1Q -23 


5 x 


1(T 15 


-S2 A131 


9 x 


1Q -24 


5 x 


10- 15 


B2 A 132 


5 x 


1Q -22 


5 x 


1(T 15 


-S2 A232 


3 x 


10 -B 


6 x 


io~ 2 


B2 A 233 


7 x 


1Q -23 


5 x 


10- 15 


-S3 A121 


5 x 


1Q -24 


4 x 


io- 16 


-S3 -^123 


7 x 


1Q -23 


4 x 


io- 16 


B3 A 132 


5 x 


1Q -22 


1 X 


10- 17 


-S3 A232 


5 x 


1Q -22 


4 x 


io- 16 


-S3 A233 


3 x 


IO -5 


3 x 


1(T 2 



Among all the BiX combinations, -B1A123, -B1A132, -S2A232 and -B3A233 are most important 
because they can provide large amplitudes via tree level diagrams. The amplitude can 
be approximated by M. ~ -BjA(tan/3sina — cos a)/ (\/2M s 2 ) , where a is the mixing angle 
between two CP-even neutral Higgs. Even though heavy sneutrino masses tend to suppress 
the amplitudes, they would relax the bounds on Bi and A more significantly, hence are 
favorable. Moreover /io should not be too small in order to make product of Bi and A to 
be below the bounds from leptonic radiative decays, i.e. | ^1 A132 1 , |-BiA* 23 |, I-B2A232I and 
I _S 3 A233 1 ~ 1-4 x 10~ 3 I /xo 1 2 [2li | . It is noteworthy that even under stringent neutrino mass 
< leV constraint, the four combinations could give branching ratios beyond IO -5 , which 
is large enough to be possibly probed at LHC. As to other BiX combinations, they can be 
from several diagrams. However, they only play minor roles and hardly give any meaningful 
branching ratio, as shown in Table 3. 

If we do not impose sub-eV neutrino mass bounds, .B1A123, -S1A132 and -B2A232 still get 
constrained by BiBi neutrino mass loops, while -B3A233 is not only constrained by BiBi 
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neutrino mass loops, but also limited by our manual "1% of B Q " bound. Comparing with 
neutrino mass < leV case, branching ratios from B 2 and B 3 increase significantly while that 
from B\ just raise a bit. 

As a matter of fact, the class of BiX combinations gives the most important contributions 
to the flavor violating Higgs decays among all RPV parameter combinations. Moreover, the 
approximation of tree level amplitudes as above could apply to h° —> ef+re and h° —> e/2+/xe 
as well. As a result, under the same parameter setting, it is expected for h° — > ef + re and 
h° — > e/2 + fie to give branching ratios with roughly the same order of magnitude as in 
h° — > /if + r/L However, it is pointed out 33] that LFV effective coupling between light 



Higgs, electron and muon could not be large because of the constraint set by two-loop Barr- 
Zee diagrams js^J on fi — > e'j. Therefore, only h° — > ef + re is expected to give branching 
ratio comparable to which of h° — > fir + r/2. 

D. Contribution from BiA^ combinations 

Under our parameterization, v4 A, s do not contribute to radiative decays such as b — > 57 



in one- loop level [2l|. Therefore, A x, s do not have known experimental constraints and, 
naively, can take any value. But the Bi parameters are limited by loop neutrino masses 
as before. Contributions from BiA x may be quite interesting since this will be like first 
experimental signature of the RPV A-parameters. However, an A x only plays its role in 
the Higgs decay through Type2No.4 diagram with a neutral scalar-charged scalar-charged 
scalar coupling. It is then expected to give a larger contribution at low charged 

scalar mass. 

Table 4. Bi A x contributions to branching ratio of h° — > [if + r/2 



RPV Parameter 
Combinations 


With Neutrino Mass 
<leV Constraint 


With Relaxed 
Neutrino Mass Constraint 


Bi A x 23 


5 x KT 11 


2 x KT 10 


Bi Ai 32 


5 x KT 11 


2 x KT 10 


B 2 A x 32 


5 x KT 11 


7 x KT 7 


B3 A x 33 


5 x KT 11 


1 x KT 7 



In our parameter setting, branching ratios from B{A X combinations can reach order of 
10 -7 at most as shown in Table 4. However, if we allow A x to be as large as 3.1TeV (or 
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higher), order of 10 6 or above is possible. Since decay rate is proportional to amplitude 
square and hence A x square, it is easy to see how branching ratio modifies as A x increases 
22\- 



E. Contribution from /x, A combinations 

All fii A combinations which can contribute to h° — > fif+rfl at one-loop level except n±\ 



123 



and /iiAi32, are constrained by their loop contributions to neutrino masses 31] . Again, single 
/ij is constrained by its contribution to tree level neutrino mass. Leptonic radiative decays 
also give upper bounds on /ij A, i.e. |/4^232| , |/4^i32| > 1/^2331 an d l^i^ml « 7.0 x 10~ 4 \fj, \ 



21] . Further bounds for single A by charged current experiments can be found in 32 1. 
Table 5. /ij A contributions to branching ratio of h° — > //f + Tp 



RPV Parameter 
Combinations 


With Neutrino Mass 
<leV Constraint 


With Relaxed 
Neutrino Mass Constraint 


Ml ^123 


5 x 10" 


-8 


5 x 10- 


-8 


^1 A132 


5 x 10" 


-8 


5 x 10" 


-8 


^232 


3 x 10- 


-12 


5 x 10- 


-8 


M2 Ax31 


2 x 10" 


-24 


2 x 10- 


-16 


M3 -^233 


1 x 10- 


-11 


5 x 10- 


-8 


M3 A 121 


1 x 10" 


-24 


5 x 10- 


-19 



Many diagrams contribute to h° — > /if + r/i process via /ii A combinations. Among which, 
TypelNo.3 and TypelNo.4 diagrams play the most important roles. The requirement of 
neutrino mass < leV still set most stringent bounds as in the case of other type combinations. 
However, //1A123 and /iiAi 32 do not give loop contribution to neutrino masses, and thus they 
are mainly bounded by the constraints from radiative leptonic decays. Generally speaking, 
large slepton masses are favorable in order to have larger branching ratios since they can 
relax the constraints from loop neutrino masses and raise the upper bounds on the A's. 

In the condition that the neutrino masses being larger than leV, the leptonic decays 
become the major sources of constraints except for the case of \ii A131 and M3^i2i- 111 an Y 
case, branching ratios from A can only achieve at most the order of 10 -8 in our analysis 
because of the stringent constraints from leptonic decays. Our results are shown in Table 5. 
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F. Contribution from the other insignificant combinations 

In addition to the above combinations, there are some other types of combination (i.e., 
Bi A', Hi fXj, Hi A', A A and A' A') which can merely give negligible contributions. Hence we 
only list the combinations which are most illustrative or give the largest branching ratios in 
each type of combinations, as shown in Table 6. Note that the types of combination which 
are not mentioned, A" A" for example, give zero contributions at one-loop level. 



Table 6. The most interesting examples in the other RPV combinations 



RPV Parameter 
Combinations 


With Neutrino Mass 
<leV Constraint 


With Relaxed 
Neutrino Mass Constraint 


^ 2 /v 333 


1 x KT 14 


2 x KT 11 


B 3 A 233 


1 x KT 15 


8 x KT 13 


H2 A*3 


2 x KT 18 


2 x KT 18 


^2 A 323 


3 x KT 18 


9 x KT 13 


A 223 


5 x KT 19 


5 x 10~ 14 


A232 A233 


2 x KT 19 


2 x KT 11 


Al21 A131 


1 x KT 15 


1 x 1(T 15 


A123 A133 


2 x 1(T 10 




A 211 A 311 


7 x 1(T 18 


9 x lO" 16 


A 222 A 322 


4 x 10~ 22 


3 x KT 15 


A 223 A 323 


4 x 10" 12 




A 233 A 333 


5 x lO" 26 


1 x lO" 15 



In the Bi A' combinations, besides the constraints mentioned before on Bi, A' also gets 
constrained by charged/neutral current experiments 32j, |35[. Bi\' combinations contribute 
to h° — > fiT + Tp, mainly via No.l and No. 2 diagrams of Type 6 and 7 (in Appendix A). 
To get better branching ratios, it is advantageous if we raise the upper bounds on Bi by 
the heavy sneutrino and neutralino masses. Heavy squark masses could also raise the upper 
bounds on A'. However, in our computation, contributions from Bi\' can not provide sizable 
branching ratios. 

As to HiHj combinations, only /i 2 /i 3 contributes to h° —> ht + tJi up to the one-loop level. 
With non-zero m , one of the neutrinos get tree level mass. However, leptonic radiative 
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decays set more stringent bound on ^2/^3 than neutrino mass does, i.e. 1/4/^3 1 ~ 3.6 x 



10- 3 lyUol 



21] . Interestingly enough, though ^2/^3 combination contributes to the decay in 



32, N|- In 



tree level, loop contribution from TypelNo.4 diagram is generally more important due to 
the smallness of neutrino masses in the loop. At any rate, ^2/^3 could only give a negligible 
branching ratio. 

Among all /ijA' 's which give non-zero contributions, some combinations are constrained 
by their loop contributions to neutrino masses 311 ] . Besides, every /ij A' is boun 
level neutrino mass constraints on /ij and experimental constraints on single A' 
this type of combinations, there is no obvious dominant diagram. Several diagrams can 
give comparable major contributions to h° — >■ \if + tJx process. Generally speaking, heavy 
gaugino masses can relax the tree level neutrino mass constraints while heavy down-squark 
masses can raise the upper bounds of A' and relax loop neutrino mass constraints, hence 
they are favorable for larger branching ratios. Unfortunately, in whole parameter space, it 
is hard for the fiiX' to give any significant branching ratios, even in relaxed neutrino mass 
constraints. 

Contributions from A A combinations are mainly from No. 3 and No. 4 diagrams of Type 2, 
6, and 7. Among all A A combinations which contribute to h° — > \if + rp,, only A232A233 and 



A121A131 are constrained by their loop contributions to neutrino masses [3 
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32 



. However, lep- 
36|. Specifically, 



tonic decays could also provide upper bounds for AA combinations 
neutrino mass constraint on A121 A131 contains a factor of electron mass, and hence is relaxed 
by the smallness of electron mass. Therefore the branching ratio from Am Am is mainly 
limited by restriction from leptonic decays. On the other hand, neutrino mass constraint 
on A232 A233 is enhanced by a r mass factor. The later gives major restriction on branching 
ratio from A232 A233 under the condition of neutrino mass < leV. The best branching ratio 
from A A we can have is of the order 10~ 10 at most. 

Many diagrams can contribute to the Higgs decay via A' A' combinations. Among which, 
TypelNo.2 and Type2No.l diagrams are most important ones. Just like the A A case, among 
all A' A' combinations which contribute to h° — > //f + t£L, only A 211 A 311 , A'222^322 an d A233A333 



as listed contribute to neutrino masses 



3l| and hence get additional constraints. Besides, 



radiative B decays and leptonic decays also give upper bounds on A' A' [32|, |36H38| . Par- 
ticularly, under the neutrino mass <leV condition, constraints on X' 211 X' 311 , A 222 A 322 and 
A233A333 are suppressed/enhanced separately by the electron, muon, and tau mass factors. 
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This makes the differences between their branching ratios. However, if we allow neutrino 
masses being larger than leV, they will constrained by leptonic decays before exceeding ex- 
perimental neutrino mass bounds. Nevertheless, A' A' type of combinations could only give 
negligible contributions to h° — >■ /if + Tp. 

IV. SUMMARY 

We have analyzed thoroughly Higgs to /if + Tp decay under the framework of minimal 
supersymmetric standard model without R parity. By means of full one-loop diagrammatic 
calculations and taking the RPV terms as the only source of lepton flavor violation, we 
showed that the branching ratio of h° — > /if + Tp could exceed 1CT 5 without contradicting 
experimental constraints. If indirect sub-eV neutrino mass bound is lifted, branching ratio 
of 10~ 2 is possible. We pull together the most interesting RPV parameter combinations and 
corresponding branching ratios in Table 7 for easy reference. Moreover, h° — > ef + re is 
expected to be able to give roughly the same order of branching ratio with that of h° — >■ 
/if + Tp from RPV terms, while h° — > ep + fie is suppressed due to stringent constraint from 
two-loop Barr-Zee diagrams. 



Table 7. Interesting contributions to branching ratio of h° — > fir + r/i 



RPV Parameter 
Combinations 


With Neutrino Mass 
<leV Constraint 


With Relaxed 
Neutrino Mass Constraint 




1 x 10" 15 


9 x KT 6 


5 3 /i 2 


1 x 1(T 13 


7 x 10~ 4 


B\ A 12 3 


1 x 1(T 5 


4 x 10" 5 


B\ A132 


3 x 10" 5 


7 x llT 5 


B2 A232 


3 x 1(T 5 


6 x KT 2 


B3 A233 


3 x 10- 5 


3 x 10- 2 


B 2 A232 


5 x KT 11 


7 x 1(T 7 


B3 ^233 


5 x KT 11 


1 x 1Q- 7 



From an experimental point of view, a typical cross section of the MSSM 125 GeV Higgs 
at 8 TeV energy is of the order lOpb. For Br(h° — >■ /if + Tp) > 1CT 5 and with a luminosity 
of the order 10 fb" 1 , this would lead to several raw fif + Tp events with almost no SM 
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background. If we allow more free parameters or a larger parameter space during our 
analysis, the branching ratios can become even larger. Together with the 14 TeV energy 
for future LHC runs, we may have more events and a better chance to probe lepton flavor 
violation, and physics beyond the Standard Model. 
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Appendix A: One-Loop Feynman Diagrams in MSSM without R Parity for Neutral 

Higgs 0° — y fir 



TypelNo.l 




TypelNo.2 




TypelNo.3 



A,i 



*3 



TypelNo.4 



T 



< 



Ui 

x 

u , 



Type2No.l 



X 



di 
X 

X 



Uk 



Type2No.2 



X 



X 
X 



T 



T 



Type2No.3 



X 

- X t 

Type2No.4 



Type3No.l 



W, 



X 

*7 



Type3No.2 




Type3No.3 




Type3No.4 




Type4No.l 




Type4No.2 




Type4No.3 




Type5No.l 
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Appendix B: Effective Couplings in MSSM without R Parity 

We list all relevant effective mass eigenstate couplings for our analysis here. Indices run 

from 1 to 10 for neutral scalars (sleptons), 1 to 8 for charged scalars (sleptons), 1 to 6 for 

squarks, 1 to 7 for neutral fermions (neutralinos) and 1 to 5 for charged fermions (charginos) 

while all dummy indices run from 1 to 3. Moreover, 

g 2 m Ui g 2 m di g 2 m { 
Vui = -JkIT. — r-: > V<k = -7^77 , and y e 



V2M W sin/3 ' ' ^2M W cos/3 ' " y/2M w cos/3 

are the diagonal quark and charged lepton Yukawa couplings, where m±s (~ m ei under 
small- Hi scenario) are mass parameters in charged fermion mass matrix 231 ] . 

Neutral Scalar (Sneutrino)-iy + -W r_ Vertices 

C = gZW + W-4L 

where 

gZ = g 2 M w (sm(3V s lm + cos(3V s 2 J . (Bl) 
Neutral Scalar-Chagred Scalar (Slepton)-^ Vertices 

where 

GZn = [(^lm ~ *^6m) ^ln ~~ (^2m — ^7m) ^2n — (^(p+2)m — *^(p+7)m) ^(p+2)n] ■ 

(B2) 

Neutral Scalar-Neutral Scalar- Z Vertices 

£ = f?S[p($)-p($)] M Z#$ 

where 

Gij = \g* (PiPli - KPh + v s 7l v s 2] - w x v 7j + v s {q+7)i v s {q+2)j - vi q+2)l v[ q+7)j ) . (B3) 

Neutral Lepton (Neutralino)-Charged Lepton (Chargino)-V^ Vertices 



Mmi 7m — o hA/ ™ 7m — o — 



)+h.c. 



where 



v/2 4m y/2 



= v u x. 



i 



--V 2i X 



Charged Lepton-Charged Lepton-Z Vertices 



C' L >y I—* + C zR i 



where 



q L = - (-1 + sin 2 ^) Ulfrj - (-± + sm 2 e w ^j U* 2l U 2j - (-± + sin 2 fl 
Cf = (1 - sin 2 ^) VijVl, + Q - sinX) V 2j V* 2i - sm 2 6 w V {q+2)j V 



i * (q+2}. 



Charged Lepton-Charged Lepton-7 Vertices 

where 



'Mo lfi 9 



$(^)*( X -) 



CJ L = U* u U l3 + C/^C7 2J + u\ q+2)i u {q+2)j 
cf = V.jVl, + v 2j v* 2i + v {q+2)3 v* (q+2)l 



Charged Lepton-Down Quark-Up Squark Vertices 



where 



u, 



dL 



1 ~ 75 c dR 1 + 75 



+ h. 



V ln i -Mn~ t ~ K CKM v 2n" Ly (j+3)m 

v* 



^nmi ~~ ~~~ u 2n u im ~i ~ u (J+2)n u hm 

fj'2 fj2 



Neutral Scalar-Quark-Quark Vertices: Down Sector 
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where 



hmi 



!hh 5 ih (v s 2m + w°j 4^(: p 



V2g 2 



V2g 2 



(k+2)m + ^(fc+7)m) 



hmi 



V2g., 



V2g 2 



(k+2)m ^(k+7)m) 



(B8) 



Neutral Scalar-Quark-Quark Vertices: Up Sector 

where 



jTruL 1-75 , rruR 1 + 75 



hrni" 



Vui 



M uR 



hmi 



V2g.. 

Vig-. 



5 ih {V{ m - iVU 
■5 ih {V[ m + W°J 



(B9) 



Charged Lepton-Up Quark-Down Squark Vertices 



where 



c uL 1 ~ 7s r „ R 1 +75 



*(tij) + h.c. 



/7uL T/ip* jt <j-\d* _|_ y^p_irip* JT <r~\d* 

^nmi * CKM 

L nmi ~~ CKM 2n pm 

g-2 



A' 

jhp-trih* jt -r-jd* 
"+" ~~ l/ CXM ty (i+2)n^(p + 3)m 

y2 



(BIO) 



Neutral Scalar-Charged Lepton-Charged Lepton Vertices 



1-75 



| ^77,77)7 



1+75 
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where 



d mn = - ^=v* 2n u ln (v s lm + w s 6m ) - -L v * m u 2n (v s 2m - w 7 j 

- y| 1/ ln C/ (i+2)n(^ , ( j +2)m ~ iV (j+7)m) ~ ^j^ V *j+2)n U U+2)n(' D 2m + iV 7m) 

+ ^*j+2)n^2n{T^lj + 2)m + ^(j+7)m) + ^(fc+2)n^(i+2)n(^ , (j+2)m + ^{j+7)m) 

Clin = ~ ^f* m V 2n {V{ m - fljgj - ^C/^V ln (P| m + iZ>?J 
1 y e . 

~ -^ U *(j+2)n V ln{^ , lj+2)m + iT> (j+7)m) ~ -^- U *j+2)n V (j+2)n(T> 2m ~ W 

+ -^^ U 2n V (j+2)n('D S u+2)m - iV S {j+7)m ) 
X* 



7m J 



+ -^^ U \i+2)n V {k+2)n{^lj+2)rn ~ tV (j+7)r 



(Bll) 



Neutral Scalar-Neutral Lepton-Neutral Lepton Vertices 



where 



' M XL 1^75 m R 1 + 75 



J v nm.r, 



6m) 



+ - (tme w x m - X 2n ) X An (V s 2m - W s 7m ) 
1 

+ - (tan$ w X m - X 2n ) X^ k+A ) n (V^ k+2 y m - iV* k+7 ) r 
Nnm n = 2^ 3n (~ ^& n @wXi n + X 2n ) {V[ m — iT>l m ) 
+ l -Xl n (tan^X* n - X* n ) (V s 2m + W 7m ) 



+ 2^-(*fc+4)n (tan^X^ - X 2n ) (D* k+2 ) m + iV° k+r}!i 



)■ 



(B12) 



Charged Scalar-Neutral Lepton-Charged Lepton Vertices 



p X L 1-75 ^ 1 + 75 



*(x°)$(0-)+h.c. 
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where 



2 



- V^X* n P< m + ^ V* m (- tan^X ln - X 2n ) V l lm 



- V2t^6 w V\ J+2)n X ln V\ q 



'i / 1 o)m ~ (^(j+2)n-^4n^ > (j+2)m - V*Q +2 y n X( j+ 4) n V 2m ) 

y2 



A, 



#2 



^*(k+2)n^{i+i)-nP\j+2)m 



Gfimn — - U* m Xl n V l 2m - C/M^(fc + 4)n^ ) (fc+2)m. + ~^^2n (tail^-X^ + X 2n ) T> 2m 

+ _ 2-^ r (fe+2)n (tan^X*^ + X 2n ) £>( fc+2 ) m 
^ e fc (tt* tV rr* v 

~ l Cy (fc+2)n^ L 4n iy (fc+5)m _ U 2n^\ 



92 

A ijk 
92 



u 



(j+2)n^(i+A)n u {k+5)m ■ 



(B13) 



Neutral Scalar-Squark-Squark Vertices: Down-Sector 



where 



d = 5A fl 1 



- l - sin 2 ^) (cos/3 ^ m - sin/3 PL) 

-V2y dq m dq Vf a V d qh V s 2m - V2y dq m dq Vf q+3)a Vf q+3)b V s 2m 
(/i* 5 pq y dq + »*\' ipq ) Vf* +3)a V d pb (V\ m + W*J 



1 



_ ^^p<Z^02%3)a^p& (^(j+2)m + ^(j+7)m) 

jpq pa (,?+3)6 ' 



(j+2)m ^(j+7)m) 



1 

71 



m dp^ipq'D d p + 3)bE > '(<l+3)a + m dqKpq^fcPpb iP(i+2)m + ^(i+7)m) 
m d P Kpq V (U3)a V (q+3)b + m d q Kpq V tb' D p l *a { V (i+2)m ~ iV {i+7)m) ■ ( B14 ) 
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Neutral Scalar-Squark-Squark Vertices: Up-Sector 

C = 9 u abm ^(u a Mu b M<j ) J 



where 



g 2 M z (2 



cos# w V 3 



frinX (cos/3 P 2 s m - sin/JDJJ ^; +3)a T^ +3)6 



+^y Uq v% M (v° m - iw r j + $ (vi j+2)m - ivi j+7)m )} v?; +3)a v; b 

+ ^Vu q V c T M [to (V2m + iW 7 J + N {Vh + 2)rn + «W>)] 
— ^^pg^"?+3)a^p& (^lm ~~ ^6m) 

--L^;p-^ +3)6 (pj m + wu . (B15) 



Cubic Neutral Scalar Vertices 

where 

4 COS (7^ 

iPlaPlb + ^IJ^lb ~ T^icPib ~ ^IcPlb ~ '^ > (q+2)aV(q+2)b ~ ^ > (q+7)aP(q+7)b) 

+ permutations in (a,b,m) . (B16) 
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Neutral Scalar-Charged Scalar-Charged Scalar Vertices 

where 

_ 1 g 2 M z . Til* Til 9?M Z I 1 . 2 Q \ QTiS Til* Til 

~ \^ir ( X " sin2 ^) [ cos/3 ^ - iV U + sin/? (V s 2m + W s 7m )\ V l * a V\ b 

Zt COS C'i^/ 

- l^TT ^ ~ Sin2 ^) [cOS/? ^™ + + sin/3 (V* m - iV'J] V[* a V 2b 

a cos v w 



'w 

' (jAfz (l-sin 2 ^) (sin/3X?i 6 + cos/3 P 26 ) (2^ +2)m + i^ +7)m ) V? q+2)a 



2 COS 6^ 

1 g 2 M z 

2cos# 



(1 - sin 2 ^,) (sin/3 Pi; + cos/3 P 2 * a ) (P^ +2)m - W s {i+7)m ) V\ q+2)b 



2cos^ C0S ^-^A - ^ ^ - ) smpV{ m V 2 * a V 2b 

+ (sin 2 w ) (cos/3 V s 2m - sin/3 V\ q+5)a V\ q+5)b 

COS t/^/ 

' Vf ^^ e9 ^ > («+ 5 ) a ^ >16 + ^7m) ~~ '^NVeq^laV \q+5)b (^2m ~~ ^7m 

+ {^pqVeq + lA^ipq) ^ > (*q+5)a'^ > \p+2)b (^lm + ^6m) 

+ {^pqUe q + HiKpq) '^ >l *p+2)aV(q+5)b (^lm _ *^6m) 

+ y| (^pqUeq + Vi^ipq) ^{q+^cPlb iPlp+2)m + ^(p+7)m) 
+ {^pqUe q + H-iKpq) ^laV(q+5)b (P{p+2)m ~ ^(p+7)m) 

+ y| m e 9 l/e 9 ^(9+2)0^26 (^+2)™ + ^( 9 +7)m) 
+ ^ m e q ye q ^2aV(q+2)b (^(^m ~ ^( 9 +7)m) 
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_ \7f mep ^ P9 ^'(9+5)a'^ ) (p+5)b (^(j+2)m + ^(j+7)m) 



-J^ m e q ^jpq^' l { q+ 2)cP{ P +2)b (P(j+2)m + ^(j+7)m) 



1 

^ \7f 777,6,3 ^P^(p+2)a^(< ? +2)6 (^(i+2)m ~ ^(i+7)m) 



^/2^pg^ > (g+5)a^'(p+2)fe (^2m + ^7m) y/2^ M "^ > ( p+2 ) a '^ , ( IJ+5 ) fe ^ 2m ^7m) 
+ 7/^^(9+5)^26 (^(p+2)m + ^(p+7)m) + 77^^*^20^(9+5)6 (^(p+2)m ~~ ^(p+7)m) 
~ 7y| yl jpg' Z: '(g+5)a :C, (p+2)6 (^(i+2)m + ^(j+7)m) 

~ 7y|^jpg' 7 ^(P+2)a'^ , (g+5)6 (^(i+2)m _ ^(j+7)m) ■ (B17) 
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